Diffusion of intrinsic localised modes by attractor hopping 
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Abstract 

Propagating intrinsic localised modes exist in the damped-driven discrete sine-Gordon chain as attractors 
of the dynamics. The equations of motion of the system are augmented with Gaussian white noise in order 
^ ' to model the effects of temperature on the system. The noise induces random transitions between attracting 

configurations corresponding to opposite signs of the propagation velocity of the mode, which leads to a diffusive 
motion of the excitation. The Heun method is used to numerically generate the stochastic time-evolution of the 

in : 

configuration. We also present a theoretical model for the diffusion which contains two parameters, a transition 
probability B and a delay time ta- The mean value and the variance of the position of the intrinsic localised 
mode, obtained from simulations, can be fitted well with the predictions of our model, and ta being used as 
parameters in the fit. After a transition period following the switching on of the noise, the variance shows a 
linear behaviour as a function of time and the mean value remains constant. An increase in the strength of the 
noise lowers the variance, leads to an increase in O, a decrease in ta and reduces the finite average distance a 
mode travels during the transition period. 



1 Introduction 

Intrinsic localised modes (ILMs) are localised excitations that can occur in a system 
due to the interplay of discreteness and nonlinearity, for instance see [1]. These modes 
exhibit internal dynamics, like oscillations or rotations of constituents of the system in 
the localisation region. Therefore, these modes also are known as discrete breathers, 
alluding to the superficial similarity between them and the breather solutions of the 
continuum sine-Gordon equation. The majority of the literature deals with strictly 
time-periodic, non-propagating ILMs. For this specific type of discrete breather there 
are rigorous proofs of existence for the conservative [2, 3] and dissipative [4] cases. From 
the proofs it is clear that the requirements for the existence of these modes are rather 
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weak, which makes ILMs quite generic, in contrast to breather solutions in the contin- 
uum. An overview of some recent results can be found in [5]. 

First predictions [6] of the existence of such modes, however, did not exclude the case 
of propagating discrete breathers and numerical results [7] soon showed the existence 
of mobile ILMs. Several aspects of such travelling modes have been investigated so far, 
for example the interaction with an impurity [8] or the effects of a bending of the chain 
along which the excitation is travelling in connection with models of DNA [9, 10]. To our 
knowledge, the diffusion of intrinsic localised modes has not been looked at yet; this is 
in stark contrast to the thoroughly investigated topic of soliton / solitary wave diffusion. 
For the latter, see for instance [11-16] and references therein. 

We study this problem in the discrete sine-Gordon chain (also known as the standard 
Frenkel-Kontorova model) under damping and harmonic driving; the equations of mo- 
tion, in dimensionless units, are: 

d^w 1 

-ttt = -7^-sm(2nu n ) + C(u n+1 + u n -i ~ 2w n ) - a—^- + Fsm(u t). (1) 
at z 2n at 

This equation describes the dynamics of a chain of identical damped and driven pendula 
in a homogeneous gravitational field, with harmonic nearest neighbour coupling. The 
quantity u n corresponds to the angle of deviation (in units of 2tt) of the n-th pendulum 
from the position of minimum energy. Certain systems of Josephson junctions can also 
be described by equation (1), see [17] for example. 

Numerical results [18] show that intrinsically localised modes exist as attractors of the 
dynamics governed by this equation. In particular, within certain ranges of the parame- 
ters C, a, u>o and F the attracting configurations are propagating ILMs with well defined 
propagation velocities, dependent on the system parameters; apart from the fact that 
the discrete sine-Gordon equation is well-known, this has been the main motivation for 
our choice of equation (1). Due to the symmetry in the equations of motion, for each 
absolute value of the velocity we can have positive and negative sign, i.e. propagation to 
the right or left, respectively. There exist regions in parameter space where more than 
one absolute value of the propagation velocity v is possible at fixed parameters. 
In order to model temperature we add Gaussian white noise terms £„(£) of strength a, 
i.e. 

(Mt)) = 

(U*)UO> =(r 2 S mn 6(t-f) (2) 

to the equations of motion (1) for the individual pendula. Such noise terms can be 
described in a clearer way by Wiener processes [19], but as the theoretical treatment of 
the diffusive motion we present in sections 3 and 4 does not make use of the formalism 
of stochastic differential equations, we do not tarry here. The relation between a and 
the temperature T is a 2 = 2aksT, where ks is Boltzmann's constant. It turns out 
that the noise causes random transitions between basins of attraction corresponding to 
opposite signs of the propagation velocity of the discrete breather, which brings about 
the diffusive motion of the excitation. 

In section 2 we present all our numerical results, including the comparison with pre- 
dictions to be derived in the next two sections. Section 3 discusses a theoretical model 



for the diffusive behaviour of the modes which uses only the transition probability 
between basins of attraction as a parameter; it is then improved in section 4 by the 
introduction of the delay time second parameter. We discuss the results and 

draw some conclusions in section 5. An appendix lists some intermediate results related 
to section 4. 



2 Numerical Results 

The simulations are done in a chain of 1000 particles with periodic boundary conditions 
(i.e. a ring). The stochastic time evolution according to (1) with the noise terms added 
is generated by the Heun algorithm (time step At = 0.01). We use parameter values 
a = 0.02, F = 0.02, lu = 0.2tt and C = 0.890 along with several values of a. The local 
energies 

En = l p * + [1 " cos ( 27ra ™)] ( 3 ) 

are used to define the position X of the localised mode as 

n=+2 

(n + n)E no+n 
X = n= ~ 2 (4) 

n=+2 v ' 

n=-2 

where the site no is chosen in the following way: The distribution of the E n sharply 
peaks around the location of the ILM. Usually, this peak has only one maximum. In 
this case, Uq is the site where this maximum occurs. Sometimes, the peak displays a 
structure maximum-minimum-maximum on three subsequent sites. In this case no is the 
site of the peak's minimum. Using (4) we calculate the position of the excitation out of 
the numerically generated configurations. Averaging over 500 realisations, we obtain the 
time evolution of the mean value of X as well as the variance of the position. In figure 1 
we show some sample trajectories. We can see that the noise causes the system to jump 
between attracting states corresponding to opposite signs of the propagation velocity 
v. At the parameter values chosen, the system is known [18] to allow for at least two 
values of |i>|, which are approximately 0.0186 and roughly twice that. As the trajectories 
shown already hint and a closer quantitative look confirms, only the former of the two 
values is involved. We cannot strictly exclude the appearance of more than just this 
one value of |i>|, but our numerical results show that if this happens, it is a rare event. 
Transitions v rts 0.0186 <-> v « —0.0186 are clearly dominant. A comparison between 
the results from simulations and the predictions of the model in section 4 is shown in 
figure 2. We observe that the mean value approaches a a-dependent value constant in 
time, i.e. after a transient there is no mean drift in the system. The variance, also after 
a transient, depends linearly on time, which is the result expected for standard diffusion 
(Brownian motion). Completely non-standard, however, is the fact that the diffusion 
constant, i.e. the slope of the variance as a function of time, decreases with increasing 
noise strength, i.e. increasing temperature. Such a temperature dependence has been 
reported for one of several diffusive processes contributing to the diffusion of kinks in a 
4 -model [13]. As we can see from figure 2, the theoretical model developed in sections 




Figure 1: Some sample trajectories for two different strenghts of the noise. Top: a 
The position of the localised mode has been calculated every 100 time units. 



= 7-10~ 5 , Bottom: a = 



12-10~ 5 . 




Figure 2: Comparison between data from simulations and results from the model of section 4. The simulation data 
are averages over 500 realisations. Left: Mean value of the position (with respect to the initial position), right: 
variance of the position. Top: a = 7 ■ 10~ 5 (fitted parameters: = 2.2 • lO -4 ,^ = 700), bottom: a = 12 • 10~ 5 
(fitted parameters: = 1.3 • 10~ 3 , t a = 280). 

3 and 4 reproduces this behaviour. The agreement for the variance is quite satisfactory, 
for the mean value it is harder to judge, because of the strong fluctuations still present 
despite the averaging over 500 realisations. 

3 A first model 

The observation of the numerically generated trajectories reveals that the effect the noise 
has on the excitation are apparently random transitions of the configuration between 
basins of attraction corresponding to opposite signs of the propagation velocity. Also, 
only one absolute value of the velocity plays a role; there is no evidence of the involvement 
of other velocity values, which in priciple are possible in the system. We consider a 
countable set of possible velocity values u» together with the probabilities Pi{t) of the 
system being in the configuration corresponding to the propagation velocity V{ at time 
t. From the obvious relation 



X(t) = J v(t')dt' (5) 
o 

for the position X(t) of the breather we obtain 

t t 



for the mean value. For the two-time correlation function C(t 2 , ti) we find 

C(f 2 ,*i) ■= (X(t 2 )*(ti)) - (X(t 2 ))(X(t 1 )) = 

*2 tl 

= v i v i J J lP( v ii 4; Vj, t[) - Pi(4)Pj(4)] dtid^ = 

M 

*2 *1 

= ViVi J J [p^(t 2 \t[) - Pi (t' 2 )] Pjit'Jdt'^t'z 



(7) 



Herein p(vi, t' 2 ; Vj, t[) is the joint probability of finding velocity Vi at time t' 2 and ve- 
locity Vj at time t[; p^jit'^t^) is the corresponding conditional probability Pi\j{t' 2 \t']) = 
p(vi, t 2 ; Vjit'^/pj^). For the case of only two velocities +v and — v the above specialises 
to 



(X(t))=v{ P + {t')dt'- / [l-p + (t')]dt' 



2 J p+(t')dt' -t 



(8) 



as for only two values p- = 1 — p + , and 

t 2 h 



C(t 2 ,t l )=v 2 




{ [p+i+it'M-pM)] P + (t[) + [p-i-it'M -p-(Q] P-(t[) 







- [p + \-{t' 2 \A) -p + {t' 2 )] P-(t[) - [p-i+teft) ~P-(Q] P+(*i)} d*id^ 2 (9) 

Up to now nothing has been said about the probabilities occuring in (8) and (9). We 
assume that there is a constant probability G per unit time for a change in the sign of 
the velocity. This leads to 

dp + 



dt 
dp- 

~dt 



e(p_ -p + ) 

Q(p+-P-) 



with the solutions 



P+(t) 
p_(t) 



i + i(p+-p_)(0)exp(-2et) 
|-|(p+-p-)(0)exp(-2e0. 



We furthermore find: 



p+i-it'M 

PM+iQt'i) 



= p.\ + (f 2 \f 1 ) = 



i i 



exp [-26 (t' 2 - t[)} 



P-\-(t' 2 \t' l ) = \ + \e^ [-20 



(10) 
(11) 

(12) 
(13) 

(14) 
(15) 



Choosing the initial condition to have positive sign of the velocity, from these expressions 
and (8), (9) we obtain 



<*(*)> = ^ [1 - ^P (-260] , £n (X(t)) = ^ 



and 



C{t 2 M) 



20 2 



exp(-20t 2 )] [cosh(26ti) - 1] 



(16) 
(17) 



From the latter expression we find for the variance of the position 

v 2 

Var [X(t)\ = 2— exp(-Gt) [sinh(0t)] 3 . (18) 

Asymptotically for long times, the variance behaves as exp(20t). This does not reflect 
the behaviour found in simulations. Thus, though due to the denominators of the prefac- 
tors in (16) and (18) an increase in the noise strength via an increase in (this relation 
between © and a is plausible and also evident from figure 1) reduces the mean value 
and partially has a suppressing effect on the variance, the results are not satisfactory 
and an improvement of the model is called for. The next section is dedicated to it. 



4 An improved model 

Individual trajectories obtained in numerics clearly show that the localised excitations 
change the sign of their propagation velocity under the influence of the noise. This 
change implies a change in the configuration of the breather, more precisely a transition 
from one basin of attraction to another. Such a transition, however, does not happen 
instantaneously. Rather, there will be a time span during which the configuration is not 
close to one of the attractors; the dynamics during this transition period may be quite 
involved and is not a subject of this paper. We attempt to roughly capture the effects of 
the non-instantaneousness of the transition by the introduction of a transition time ta, 
during which the velocity of the breather is 0, and during which no further transitions 
may be initiated. This means that if a transition from the +v state is initiated at time 
t = 0, then the velocity of the localised mode immediately acquires the value 0, which it 
will retain up to t = ta, when the velocity jumps to —v. In particular, we do not take 
into consideration 'failed' transitions, i.e. jumps of the velocity from +v to and back to 
+v. It is unclear whether taking into account such processes would improve the model; 
reality is much more complex and can involve all kinds of trajectories in velocity space. 
It certainly would increase the number of parameters and assumptions in the model; 
therefore we confine ourselves to the simple model where each initiated transition after 
time ta ends in the attractor corresponding to the opposite sign of the velocity. We keep 
the assumption of the previous section that there is a constant probabibility per unit 
time for a transition to be initiated. 

The introduction of the 'delay time' ta limits the maximum number m of jumps occuring 
in time t to m — \t / ta\ + 1, where [x] denotes the integer part of x. If we consider a 
trajectory with n jumps up to time t (note that this implies t > (n — 1)ta), the jumps 
occuring at times Ti, T 2 , . . . , T n , with < T\ < T 2 — t a < ■ ■ ■ < T n — t a , we obtain for 
the distance traversed by the ILM: 

X(t) =v[T 1 - (T 2 - T x - T A ) + (T 3 -T 2 -t a )-... 
Tv(n)(T n - T n _x - t a ) + n(n)(t -T n - t a )\ = 



- 2 J2^k)T k + Hn)t + ^-^-T A 
k=i 



(19) 



where 7r(n) is the parity of n, i.e. ir(n) — +1 if n is even and 7i(n) — — 1 if n is odd. 
The expression (19) holds if T n < t — r A \ in the case T n > t — r A we find 



X(t) = v 



2j2^k)T k -n(n)T n + 1 -±^- 



ta 



k=i 



(20) 



because the excitation does not move after the n-th jump. The temporal probability 
density of n jumps occuring at the times T±, . . . , T n is in the case T n <t — ta- 

p(T n , . . . ,Ti) = exp(-6T 1 )0 exp[-0(T 2 - T x - r A )\ . . . 
. . . exp[-0(T„ - T n _! - TA )]e exp[-G(t -T n - r A )\ = ©" exp[-6(t - nr A )\ (21) 

and analogously in the case T n > t — t a : 

p(T n , ...,T 1 ) = 0™exp[-0(T„ - (n - l)r A )}. (22) 

Using the equations (19-22) the mean value of X(t) can be calculated as 



m— 1 



(X(t)) = S (t) + [*VW + S ^(t)) + S 3 ,m(t) 



71=1 



with 



S (t) = vtexp[-et] 
corresponding to no jumps, and, using the step function H, 

t—TA T n -T A 

S hn (t) = H(t-nr A )ve n exp[-&(t-nr A )} J dT n J dT, 



n— 1 



(n-l)r A (n-2)r A 



Tz—TA Tl—TA 

... / dT 2 / dT, 



u{k)T k + 7r(n)* + -^r A 



k=i 



(23) 



(24) 



(25) 



Tn-TA 



S 2 , n (t) = H(t - nr A )vQ n J dT n exp[-0(T n - (n - l)r A )\ J dT n _ x . . . 

t—TA {n-2)T A 

Tz-ta T 2 -ta 



dT 2 J dTx 

TA 



n-1 



E,, s „ , Nm 1 + 7r(n) 
7r(A;)T fc - 7r(n)T„ + -^r A 



k=l 



and 



(26) 



x 



S 3 , m = if (i - (m - l)r A ) [1 - if (i - mr A )] vQ m x 

« T m -T A Tz-ta T 2 -t a 

J dT m exp[-Q(T m -(m-l)r A )] J dT m _! . . . ^ dT 2 ^ dT\ x 

(m-l)rA (m-2)rA 

m— 1 







X 



-2^7r(A;)T fc -7r(m)r m + i±^ 



<~.4 



k=l 



• (27) 



Several expressions are useful in evaluating these quantities as well as others to be 
introduced below. They are gathered in an appendix. After some algebra we arrive at 

,[6(f - nr A )] n+1 



S 1 , n (t) = H(t - nr A ) l + ' n 2 {n) ^ ) exp[-6(t - nr A )Y + 



(28) 



1 — 7V\Tl) V 

S 2 , n = H(t - nr A ) — exp[-9(t - nr A )} x 

E T^yy {m - nr A )r i - m ~( n - l^re^} 



(29) 



S 3 , m = H(t - (m - l)r A ) [1 - H(t - mr A )\ 1 ! M X x 



x\l-exp[-Q(t-(m-l)r A )]J2 



2 

m—l 



1=0 



[Q(t - (m - 1)t a )] 
(m-l)\ 



(30) 



The expressions for S^n an d £3,™ are complicated because of the polynomials appearing; 
in all the cases Si >n ,S 2 , n and S 3)in , however, we clearly see the 'retardation' caused by 
the introduction of t a . 

In a similar way, the expectation of X(t) 2 can be written as 

m— 1 

{x(tf) = Q (t) + £ [Q 1>n (*) + g 2 , n (0] + Q 3> m(0 (31) 



with 



n=l 



Q (t) = w 2 t 2 exp(-9t) 

[6(f - nr A )] n+2 



Qi, n (t) = H(t - nr A )—exp [-Q(t - riT A )} 



x 



(n + 2)! 



x 



^(n + 2) + -^(n + 1) 



(32) 



(33) 



n+l 



g 2l n(*) = #(* - nr A )^ exp [-6(t - nr A )\ { - ^)] n 



+i-fe 



fc=0 



X 



1 + vr(n) 



_[0 (t _ (n _l )rA) ]«+^e-e^| x 

1 — ir(n) n + l 
2 (n+l - k)\ 



n 



2 (n + l - fc)! 



+ 



(34) 



Q 3 , m (t) = H(t - (m - l)r A ) [1 - if (* - mr A )] — x 

1 + ■Kim) 1 - ■Kim) . 1 . 
2 m + 2 (m + } 

£ [9( t _ ( m _ l)r A )] m+1 - fc exp [-6(f - (m - l)r A )\ x 

k=0 



x |^l-exp[-e(t-(m-l)r A )]) ( 



x 



1 + 7t(to) 



(m + 1 - k)\ 



+ 



1 — n(m) m+1 
~2 (m + 1 -k)\ 



(35) 



From these results we can of course obtain \/ar[X(t)] = (X(t) 2 ) — (X(t)) 2 , but as the 
expressions shown above are already rather involved, we do not carry out this step 
explicitly. However, from the above results the mean and the variance can easily be 
plotted and compared with the results from simulations, as done in figure 2. The results 
from fits for values 5- 1CT 5 < a < 15- 1CT 5 clearly show that ta decreases with increasing 
noise strength (roughly as ta oc i/(/cbT)), whereas increases strongly, approximately 
as oc exp[— /i/(/c B T)] with some constant /i. Due to the strong fluctuations, clearer 
results are not possible; therefore we also do not show graphs. 

The behaviour for short times t < Ta, Qt <C 1, can be obtained from S 3j i and Q^\. We 
find 



(*(*)> = £[l-exp(-et)]«£ 



et - \ m 2 



(36) 



and 



Var[X(t)\ = ^ [1 - 20t exp (St) - exp (-26*)] » ^ (0t) 3 . (37) 

A completely different approach, which we will therefore present in a separate publica- 
tion, easily yields the long-time (t 3> 1/0) behaviour. The results are (X) = t>/(20), 
which agrees with the result from section 3 independently of ta and Var[X(i)] = Dt + c 
with 

v 2 v 2 

D = 7 TTi c ~ 7^7- (38) 

<3> 2 (TA + ±) , 20 2 1 ; 

Here, even in the case ta — > 0, an agreement with the simpler model of the previous 
section, where the variance shows exponential behaviour, is not achieved. 



5 Discussion and Outlook 

As we have seen, the diffusive behaviour of an intrinsic localised mode in the damped- 
driven discrete sine-Gordon chain, when Gaussian white noise is coupled to the system 
to model the effects of temperature, is due to random transitions between attracting 
configurations of the system that correspond to opposite signs of the propagation veloc- 
ity of the mode. In this paper we have developed a model that satisfactorily describes 
this diffusion process. At a preliminary stage, the model only uses the transition prob- 
ability per unit time, 0, as parameter. This approach is not sufficient; it is improved 
by the introduction of a second parameter, the delay time ta- Both parameters are de- 
termined from fits to the results from simulations. These fits show that with increasing 
temperature ta decreases and increases. The long time behaviour for both models 
is the same for the expectation of the position of the intrinsic localised mode, but not 
for the variance of the position, not even in the case of vanishing delay time. This is 
due to a fundamental difference between the two approaches we have presented. In the 
one-parameter model we have only taken into account the probabilities of finding the 
configurations in one of the attracting states or the other. The two-parameter model, on 
the other hand, follows single trajectories through time and considers the probabilities 
for repeated jumps between the basins of attraction to occur in one and the same tra- 
jectory. The nonvanishing delay time, during which propagation of the excitation and 
further jumps are forbidden, also introduces a non-Markovian element in the diffusion 
process. The model does not make use of any specific characteristics of the system; it 



only requires the existence of the attracting states and the possibility to excite transi- 
tions between these states by noise. The details of such a transition will, we suspect, be 
sensitive to the particular system, but in the model that we have presented in section 4 
the details are 'hidden' in the delay time t a . Thus the origin of deviations between the 
prediction of the model and the result of simulations, which show in the transient period 
between the very early stages of the stochastic evolution and the long time diffusive 
regime, lies in the dynamics of the transition from one attractor to the other. 
A first step in the development of a more detailed model would be a calculation of the 
parameters and ta from the system parameters C, a, F, ujq and the noise strength a. 
A further step would be to go beyond the simple picture of 'delayed jumps' that is at 
the heart of our model and to include the dynamics in full. The latter step then should 
also be able to describe the behaviour in the transient time regime. 
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Appendix 

We list here several expressions that arise in the evaluation of the mean value of X and 
of X 2 . M in the equations below is not to be confused with the maximum number of 
jumps m. 

Tm-ta t m _ x -t a t 3 -ta t 2 -t A Af _ l 

/ dZi,., / dT M _ 2 ... / dT 2 / dT. ^-g:^' " (39) 

(M-2)r A (M-3)t a t a 

For 1 < k < M - 1 

T M ~T A T M _ X -T A T 3 -T A T 2 ~T A 

J dT M -i J dT M _ 2 ... J dT 2 J dT!T fc = 

{M-2)t a (M-3)t a r A 

k [T M - (M - l)r A f [T M - (M - l)^]"- 1 

k M\ +[k 1}TA (M-l)! (40) 



Tm-ta Tm-i—ta Tz-t a T 2 -t a 

dT M _ 2 ... J dT 2 j dTiTfc 

TA 



M-TA J-M-l- 

I / 



(M-2)t a (M-3)t a 

k(k + DEkzigz^ + 2r Mk - 1} PW^zifclL + 



(M+ 1)! 



(k-1 



.2_2 



Ml 

[T m -(M-\)t a 
(M-l)! 



M-l 



(41) 



For 1 < k < I < M - 1 

Tm-ta Tm-i—ta T 3 —ta T 2 -t a 

J dT M _! J dT M _ 2 ... J dT 2 J mT h Ti = 

(M-2)t a (M-3)t a t a 

(H-i-^i ^-g:^ 

E l-7r(n) 
*(*0 = g 

fe=i 

fc=i ^ 7 

m/ 1 + tt(w) n 1 - 7r(w) / n + 

E ^ = 2 + 2~ ) - 

i=k+i v 7 

l + 7r{k)k l-n(k) ( fc + 1 



fc=i i=k+i v 7 v 

n-l n 1 1 , / \ 1 

fc=l «=fc+l 

E/,m V""^ /,x 1— 7r(n) / n 2 1\ l+7r(n). n 2 

fe=i i=k+i " v 7 



n— 1 n 



£ tt(/)A;/ = 

fc=l Z=fc+1 



,3 



6 



n 2 


n ^ 






12 J 


+ 


8 


n 2 


n 












6 + 
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